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The Rocketdyne Multifunction Tester is a general purpose test apparatus which utilizes 
axial and radial magnetic bearings as shaft excitation devices. The tester is modular in de- 
sign so that different seal and bearing packages can be tested on the same test stand. The 
tester will be used for rotordynamic coefficient extraction, as well as life and fluid/material 
compatibility evaluations. Use of a magnetic bearing as a shaft excitation device opens up 
many possibilities for shaft excitation and rotordynamic coefficient extraction. In addition 
to describing the basic apparatus, this paper will discuss some of the excitation and extrac- 
tion methods considered, and detail the chosen method. Some of the excitation methods 
to be discussed include random, aperiodic, harmonic, impulse and chirp. 

INTRODUCTION 

Recent trends in the design of high-performance turbomachinery for advanced pump-fed 
liquid rocket engines have moved away from rolling element bearings in favor of fluid film 
elements. The potential payoffs are longer life, higher admissible DN values, and more con- 
trol over turbopump rotor dynamics. Hydrostatic bearings and annular seals appear most 
promising. In effort to support the new designs, Rocketdyne has designed, and is currently 
fabricating, a multifunction bearing and seal tester. The goals are to measure accurate 
rotordynamic coefficients, conduct durability tests, and evaluate different materials, all in 
an environment simulating real cryogenic turbopump operating conditions. Measurement 
of rotordynamic coefficients is accomplished via the use of a radial magnetic bearing as 
a shaft motion exciter. Test bearing (or seal) reaction load is measured in the support 
structure. The versatility of the magnetic bearing as an exciter is the key feature of this 
tester. 

The purpose of this article is to describe the basic methods which will be used to excite 
the shaft, acquire the data, and reduce the data. All three of these tasks are to be 
performed by a single computer. The tester is designed to run the test bearing (or seal) at 
a constant operating point, and use the magnetic bearing to excite the shaft in a manner 
suitable to permit identification of the bearing’s (or seal’s) complex impedance functions. 

‘This work was partially supported by United States Air Force contract F04611- 
86-C-0103. 
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The impedance functions then yield the desired rotordynamic stiffness, damping and mass 
coefficients. 

The application of a magnetic bearing in this manner is novel, thus there is a degree of 
uncertainty as to how well it will function, and what the optimum input signal will be 
(i.e., random, aperiodic, harmonic, swept sine, impulse, step, chirp, etc.). The versatility 
of having a completely arbitrary waveform generator is, therefore, most important. Also, 
it is desirable that the input signal delivered to the magnetic bearing be adaptive so that 
optimum signal to noise ratios can be achieved. 

The key constraint in this test program is length of test time. Many of the tests are 
limited to about 5 minutes each by virtue of the blow down facility being used. In all, 
approximately 35 such tests are currently planned. These 35 tests are divided among 
different seals, different bearings, and different working fluids. Each 5 minute test must 
encompass a multitude of operating points (combinations of shaft speed, pressure drop, 
inlet swirls, eccentricities, etc.). Presently, only about 20 seconds of test time are being 
allotted to each operating point for measurement of rotordynamic coefficients. The chosen 
input signal must make optimum use of this limited amount of time. 

DESCRIPTION OF APPARATUS 

The multifunction tester, configured with an an n ular seal, is shown in figure 1. The tester 
consists of the basic tester and a module containing the test article. The basic tester rotor 
is supported on one end by a duplex ball bearing and on the other by the magnetic bearing. 
The module contains the test article only. When assembled, the main rotor supports are 
the duplex ball bearing and the test article, use of the magnetic bearing as a support is 
optional as discussed by Hawkins, et al. (ref. 1). 

The main function of the magnetic bearing is to generate asynchronous motion of the rotor 
at the test article. The reaction force at the test article is measured using a calibrated 
cylindrical spring called a flexmount. The deflection of the flexmount under load is mea- 
sured using four Bently REBAM proximity probes. The acceleration of the flexmount is 
measured in two orthogonal directions using PCB accelerometers. The motion of the rotor 
at the test article is measured using Bently proximity probes. 

The initial test plan calls for testing of a 152.4 mm (6.0 inch) diameter annular seal, a 
88.9 mm (3.5 inch) diameter hydrostatic bearing, and a 57 mm (2.24 inch) ball bearing in 
LN2 and actual cryogenic propellants. The annular seal will be tested at three different 
values of inlet swirl. The tester will operate at speeds up to 24,000 cpm and supply 
pressures to 172.4 Bar (2500 psi), for this test series. The maximum tester design speed is 
30,000 cpm. The geometries for the test articles will reflect those envisioned for production 
rocket engine turbomachinery. This requires radial clearances in the range of 0.076-0.127 
mm (0.003-0.005 inches), whereas all test results to date are for radial clearances in the 
range 0.254-0.508 mm (0.010-0.020 inches) and at lower shaft speeds. 
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FUNDAMENTALS OF ROTORDYNAMIC 
COEFFICIENT MEASUREMENT 

The standard relationship between bearing displacements and forces in the time domain 
is expressed as follows: 


[M\q+[C\q+\K]q = F (1) 


Rotors operating at steady state execute motions which are strongly periodic (i.e, har- 
monic). This lends itself to a frequency domain representation. Also, from a test point of 
view, it is more convenient to work in the frequency domain. Thus, transforming equation 
(1) to the frequency domain and employing complex notation (equivalent to amplitude 
and phase) results in 


(K„, + iC„w - 
(ir„„ + iC yx u - M, x ui 2 ) 


(K„ + iC„ w - M„X) 1 J 1 \ = J M 
(Ky, + iC„U - MyyW 2 ) J \ B / \ Fy J 




(2a) 


(26) 


In the usual sense, the physical displacement in the x direction is the real part of the com- 
plex expression xe iut and similarly, for y, F x and F y . The complex impedances, Z, come 
directly from measurements of x, y, F x and F y , all presumably as functions of frequency. 
The displacements x and y are usually obtained easily with non-contacting displacement 
probes. The forces F x and F y are more difficult, and may be measured directly with spe- 
cial load cells, or inferred indirectly from some other known quantities such as pressure 
distribution, a hammer impulse, intentional unbalance, or the snap of a suddenly released 
static load. Once the displacements and forces are obtained as functions of frequency, 
the exact procedure for obtaining the impedances Z, and in turn, the rotor dynamic co- 
efficients, depends on how the displacements and forces were generated (i.e., the overall 
test approach). The reduction procedures for several test approaches considered for the 
multifunction tester will now be described. 


Single Circular Orbits - A very natural approach is to use the magnetic bearing to 
impart a circular whirl orbit to the rotor, and vary its frequency over the frequency range 
of interest (rotor speed held constant). In this case y = xe~'* (counterclockwise orbit) at 
each value of frequency giving the following: 


Z mm + Z xy e~^ = ~ = K XX + C xy u - M xx w 2 - i{K xy - C xx u - M xy u 2 ) 
% 

F 

Z yx + Z„e-‘* = ~ = K yx + C yy u - M yx u 2 - i(K yy - C yx u - M yy u 2 ) 
x 


( 3 ) 


The x, F x and F y 
pletely determine all 4 


data obtained for a single circular orbit are not sufficient to corn- 
complex Z’s. However, the special shape of the orbit does permit 
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their partial determination by grouping of their real and imaginary components as shown. 
The 4 components (real and imaginary parts of F x /x and F y /x) can be plotted versus 
frequency, and curve fits yield the 12 desired rotordynamic coefficients. Note that all 12 
coefficients must be real which precludes the detection of mechanisms like hysteretic damp- 
ing. Note also that there is no benefit to assuming skew symmetric stiffness, damping and 
mass matrices. The only information obtained as to the adequacy of the assumed KCM 
model is via digression from a pure quadratic. This test approach is used by Adams, et al. 
(ref. 2). 

Dual Circular Orbits - The magnetic bearing is capable of producing either forward 
or backward circular whirl orbits at each frequency (not simultaneously). This forms 2 
independent load cases which, when used together, completely determines all 4 complex 
Z’s at each test frequency. 

'Xj j// 0 0 

0 0 x f 27/ 

x b y b 0 0 

.0 0 x b y b 

The 8 components are plotted versus frequency, and curve fits yield the desired rotordy- 
namic coefficients which can now all be complex. This test approach was used by Kanki, et 
al. (ref. 3). Note that it is better to combine the forward and backward data to completely 
determine the Z’s over the frequency range from f iow to f high , rather than keeping them 
separate and using equation (3) to partially determine the Z’s from — / high to +f h i g h (see 
Ohashi, et al. (ref. 4) and Jery, et al. (ref. 5)). This same conclusion was also reached 
by Bolleter, et al. (ref. 6). Note also that forward and backward circular orbits taking 
place at the same time form an ellipse, but a single ellipse is a single load case and does 
not completely determine the Z’s (see below). 

The technique of actually plotting and inspecting the impedance functions versus fre- 
quency, followed by a curve fit, provides insight into the mechanisms at work, and also 
experimental scatter. The importance of this visualization apparently has been overlooked 
by some researchers. As opposed to plotting and curve fitting the impedance functions, 
the force and displacement data for some or all frequencies can be assembled into a single 
matrix equation as follows (i.e., rearrangement of the terms of equation (2a)): 


Z xx 

%xy 

Z„ x 


J yy 



f = forward 
b = backward 


(4) 


[displacement data matrix]{KCM vector} — {force data} (5) 

[4 * # frequencies by 12]{12 by 1} = {4 * ff frequencies by 1} 

The 12 real rotordynamic coefficients are computed directly by either using a subset of 
the data to write exactly 12 equations (Adams, et al. (ref. 2)), or using all the data and 
performing a least squares solution (Nordman, et al. (ref. 7), Goodwin, et al. (ref. 8), 
and Burrows, et al. (ref. 9)). The latter approach will produce rotordynamic coefficients 
identical to those produced by plotting and curve fitting, but without the benefit of seeing 
this graphically. 
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Another approach to bypassing the impedance functions and proceeding directly to 12 
unknown rotordynamic coefficients is to use a recursive estimation scheme on time domain 
data (Stanway (ref. 10), Sahinkaya, et al. (ref. 11), and Ellis, et al. (ref. 12)). This should 
yield a result similar, if not identical, to the frequency domain approach of equation (5). 
Among its advantages are transformation to the frequency domain with its associated 
windowing functions is not needed, rotordynamic coefficients can conceivably be computed 
on-line during test (i.e., like a filter), and it is an effective way to deal with time varying 
parameters (Goodwin and Payne (ref. 13)). Again, however, there is will be no graphical 
check on the validity of the KCM model. 

Single Elliptic Orbit - The magnetic bearing can also be used to impart elliptic motion 
to the rotor, and vary its frequency over the frequency range of interest (rotor speed held 
constant). Data for a single elliptic orbit contains the same amount of information as a 
single circular orbit, but must be utilized differently because of its different form. For 
example, to use the technique of equation (3), y/x must be constant. The displacement 
orbit must have a constant ratio of major to minor axis length, and must have constant 
orientation. It is probably not practical to satisfy these conditions. Instead, the approach 
of equation (5) could be used to directly compute 12 real rotordynamic coefficients. In 
this case it is desirable that the displacement ellipse be highly variable as a function of 
frequency to get the best definition of the rotordynamic coefficients. In fact, the further 
removed from a circle the data becomes, the more important it is that the ellipses be 
highly varied. It can be shown that the limiting case of a constant, purely translational 
orbit would render the matrix of equation (5) singular. When solving the least squares 
form of equation (5) it is best to use a Singular Value Decomposition (SVD) algorithm 
(Forsythe, et al. (ref. 14)) to get a direct indication of the quality of the solution (i.e., 
the matrix condition number). One can also use the statistical parameters employed by 
Burrows, et al. (ref. 9) to help quantify the accuracy of the least squares solution. It is 
interesting to note that some simple numerical experiments on equation (5) have shown 
that varied elliptic data can easily yield a solution twice as good as circular data. In 
exchange for the more precise solution, however, it is no longer possible to plot impedance 
functions versus frequency to provide a visual check of the validity of the KCM model. 

Another approach to utilizing data for one ellipse per frequency is to assume the im- 
pedances to be skew symmetric (i.e., Z xx = Z yy and Z xy = —Z yx ). In this case equation 
(2b) becomes (see Bolleter, et al. (ref. 6)): 



As long as a is not ±i (i.e., not a circle), the two complex impedances can be computed 
using Cramer’s rule, and plotted versus frequency. The plots provide a check on the KCM 
model, and curve fits yield 6 complex rotordynamic coefficients ( K xx , C xx , M xx and K xy , 
C xy , M xy ). The coefficients would be expected to be real, but need not be. For example 
Bolleter, et al. (ref. 6) detected hysteretic damping (complex K xx ) in boiler feed pump 
impellers using this approach. Note that the orbit must be elliptic. As a approaches ±t 
the solution for Z xx and Z xy becomes ill-defined, and is undefined when o = ±*. The 
following expression is an amplification factor for uncertainties in both x and y: 
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2(|a| - 1) 


( 7 ) 


If there is 5% uncertainty in x then there is f*5% uncertainty in Z xx and Z xy . The need 
to avoid circles in this case is apparent. One other potential drawback of assuming skew 
symmetric impedances is what happens when this assumption is violated. For example, if 
K xx equals 90% of K yy this could manifest itself as a bias in K xy , or produce seemingly 
significant amounts of C xy when in fact there is none. 


Dual Elliptic Orbits- The magnetic bearing could be used to produce two different 
elliptic orbits at each test frequency (not at the same time) . This forms two independent 
load cases, and is analogous to the case of forward and backward circles described earlier. 
In fact, the forward and backward circles is a special case of 2 ellipses. When using two 
ellipses it is imperative that they be independent so that the 4 Z’s can be computed 
employing equation (4) for each frequency. The condition number (Forsythe, et al. (ref. 
14)) of the matrix in equation (4) quantifies the degree of independence. The following 
expression due to Zhang, et al. (ref. 15) can be used to calculate the condition number at 
each frequency: 


COND = + (8) 

B - \JB 2 - 4\x 2 y ± - y 2 x 1 \ 

where B = \x x | + \x 2 1 + \y t \ + |y 2 1- The condition number is similar to / of equation (T), 
and in this case is an upper limit on amplification of error. The minimum, or optimum, 
value is 1 which occurs only with orthogonal translational orbits of equal amplitudes. The 
maximum value is infinite which occurs with identical orbits (can have different sizes). A 
pair of forward and backward circular orbits of equal radius yields a condition number of 
5.83. Adequate solutions to equation (4) can be plotted versus frequency for subsequent 
inspection and curve fitting. Excessive condition numbers will make the impedance curves 
appear ’’noisy” regardless of the actual signal/noise ratios. This latter fact is easily over- 
looked. When the condition number is too high one must resort to the techniques for single 
elliptic data. Numerical experiments presented by Zhang, et al. (ref. 15) dictate that a 
practical cutoff on condition number may be about 20 to 30. 

It has been suggested that the 4 stiffness and 4 damping values (no inertia) can be obtained 
solely from unbalance response data at a single speed (Lund (ref. 16), Zhang, et al. 
(ref. 15), Lee, et al. (ref. 17) and related paper by Verhoeven (ref. 18)). This fits the 
two ellipse case for a single frequency. Either two identical bearings on one shaft must 
experience different elliptic orbits simultaneously under the action of a single unbalance 
distribution, or one bearing must experience different ellipses at the same speed under 
two different unbalance distributions. Either case can apparently be achieved in practice 
with a carefully designed experiment. The degree of success will depend not only on the 
accuracy of the bearing displacements and forces, but also on the value of the condition 
number of equation (8). 
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CURRENT TEST APPROACH 

In light of the preceding discussion, it has been decided that, if possible, a test approach 
would be used which enables plotting of impedance functions versus frequency without 
assuming skew symmetry. This implies either the single circular orbit approach, or dual 
elliptic orbit approach. The dual elliptic approach is preferable as it completely deter- 
mines the impedance functions, but may require twice the test time as the single circular 
approach. As mentioned earlier the length of test time is severely limited. So arbitrarily 
large numbers of samples will not be available for averaging out system noise. As a first 
step we will attempt to achieve purely translational rotor motion first in one axis and 
then the other. This would provide the most accurate impedance determination possible 
according to equation (8) . The system may, however, prove to be too noisy such that there 
is not ample time to get enough averages when executing the two load cases sequentially. 
Burrows, et al. (ref. 9) and Yasuda, et al. (ref. 19) applied their load cases simultaneously 
via alternating harmonics and statistical independence, respectively. If necessary, we will 
try the method of Burrows, but at this time it is unclear whether these are viable ways 
to overcome a problem of insufficient averaging time. If the impedance measurements still 
suffer from excessive noise, we will resort to the single circular orbit approach. Should it 
happen that the magnetic bearing cannot satisfactorily enforce a circular orbit, we shall 
use the method of single elliptic data. Depending on the overall success, or lack thereof, 
of the current test program, future testing may employ recursive time domain methods 
mentioned earlier. 

EXCITATION SIGNAL 

The radial magnetic bearing serves as a non-contacting, dual-axis electrodynamic shaker. 
Its power supply has an input jack for each axis such that it can perform as 2 independent 
shakers. In actuality the magnetic bearing system is unproven in this regard, and at this 
time it still remains to be seen how well it will perform this function. Our goal is to 
measure the transfer functions of bearing force to bearing displacement. Since test time 
is at a premium, the input signal should provide a good combination of accuracy and 
rapidity. Another consideration is the frequency range of interest will span up to 4 system 
resonances which will be functions of pressure drop, shaft rpm, etc. The magnetic bearing 
system itself also has highly nonlinear frequency response characteristics (Hawkins, et 
al. (ref. l)). It is desirable that the input signal can be spectrally shaped in an adaptive 
fashion so a flat displacement spectrum can be achieved on output. This will enable higher 
signal to noise ratios over the entire frequency range without overloading the system at 
the resonances. A flat displacement spectrum also makes it easier to investigate amplitude 
dependent nonlinearities by driving the system to progressively higher response levels. 

Single frequency steady state sinusoidal testing is unquestionably the best way to fit a 
linear model to a real system (Herlufsen (ref. 20)), but the long test time precludes its 
use here. Fortunately, work by Bolleter, et al. (ref. 6), Yasuda, et al. (ref. 19) and 
Chang, et al. (ref. 21) indicate that a multifrequency input signal gives virtually identical 
results as single frequency testing. Yasuda’s testing was steady state as they used random 
input signals. Bolleter’s testing was transient as he used a rapid sine sweep, or chirp. 
Chang’s testing was also transient as he used impact. The article by Herlufsen (ref. 20) 


369 



discusses some of the aspects, and different types, of steady state and transient testing. 
Since we intend to modify the input signal during test in an adaptive fashion to achieve a 
flat displacement spectrum, a transient signal is the logical choice. 

Step (Morton (ref. 22)), impulse (Nordman, et al. (ref. 7)) and chirp (Bolleter, et 
al. (ref. 6)) are examples of transient input signals. Bolleter obtained excellent results 
using a chirp input. Also, the spectrum of a chirp can be easily shaped to virtually any 
desired form. Papers by White, et al. (ref. 23) and White (ref. 24) provide the technical 
aspects of using a chirp input signal. The chirp is thus our chosen input signal (see figure 
2). Our test approach will be to input a series of chirps along one axis including the 
modification of each successive chirp to converge on a flat displacement spectrum. This is 
then repeated along the other axis. Very elongated elliptic orbits should result as long as 
the rotordynamic cross- coupling is not too strong. This data will then be used with the 
dual elliptic reduction method described earlier. 

DATA ACQUISITION AND PROCESSING 

A digital computer will be used to sample, and store on-line, all data channels necessary to 
produce rotordynamic coefficients. Prior to testing, a common broadband random signal 
will be input simultaneously to all channels of the data acquisition system. Interchan- 
nel transfer functions thus obtained will be used as correction factors to match the gain 
and phase of all channels, including filters. The same computer system will be used to 
compute and output the chirp signals. The computer system is capable of performing an 
on-line spectral analysis of the displacement response, and modifying each successive chirp 
accordingly. 

The complex quantities x, y, F x and F y are needed for equation (4) so that the impedances 
can be calculated. The effects of system noise can be more effectively reduced if one works 
wi fV tranr r '" funct'-mr computed from spectral densities (see Bendat, et al. (ref. 25) or 
Halvorsen, et al. (ref. 26)). Thus in practice the computer will be used to compute cross 
and auto spectral densities from each series of chirps. These are used to compute transfer 
functions F x /x, F x /y, F y /x and F y jy as functions of frequency. For each frequency these 
8 complex values (4 from the x chirps and 4 from the y chirps) are used with equation (4) 
to compute the complex impedances, Z. All calculations are to be done post-test by the 
same computer used for acquisition. Hard copy plots of the impedance functions, along 
with curve fits for rotordynamic coefficients, are the end product of each test. 

SUMMARY 

A multifunction test apparatus has been described which will be used to measure rotor- 
dynamic coefficients of bearings and seals. A radial magnetic bearing will be used to load 
the shaft statically and excite dynamic motion of the rotor in a controlled fashion which 
will enable complete definition of the test element’s rotordynamic impedances as functions 
of frequency. Transient testing will be performed with chirp input signals being applied 
separately to the inertial x and y axes of the rotor. The chirp input signals will be modified 
on-line to achieve flat displacement spectrums at the test element. The test approach and 
the data reduction methodology were described in detail. The testing environment closely 
matches real rocket engine turbopump operating conditions in speed, pressure, clearance, 
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temperature, working fluid, etc. Overall, this is considered a very ambitious test program, 
technically. Initial testing is currently expected to take place in the su mm er of 1990. 

NOMENCLATURE 


a 

y/x 

B 

|*i | + |*a I + |yi | + life | 

COND 

matrix condition number 

F 

{F.,F,} T 

F*,Fy 

real test element force 

F F 

x x 5 x y 

complex test element force 

flow ) f high 

limits of test frequency range 

f 

error amplification factor 

K 

test element stiffness matrix 

C 

test element damping matrix 

M 

test element inertia matrix 

KCM 

stiffness, damping and inertia 

q 

{x,y} T 

x,y 

real displacements of test element 

x,y 

complex displacements of test element 

Z 

complex impedance of test element 

i 

\/=l 

u 

frequency of harmonic motion (rad/sec) 

|x| 

magnitude of complex number x 
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FIGURE 1) Schematic of multifunction tester. 
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